Electronic Structure of Pyrochlore Iridates: Prom Topological Dirac Metal to Mott 
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In 5(i transition metal oxides such as the iridates, novel properties arise from the interplay of 
electron correlations and spin-orbit interactions. We investigate the electronic structure of the 
pyrochlore iridates, (such as ¥2^207) using density functional theory, LDA+U method, and effective 
low energy models. A remarkably rich phase diagram emerges on tuning the correlation strength 
U . The Ir magnetic moment are always found to be non-coUinearly ordered. However, the ground 
state changes from a magnetic metal at weak U , to a Mott insulator at large U . Most interestingly, 
the intermediate U regime is found to be a Dirac semi-metal, with vanishing density of states at 
the Fermi energy. It also exhibits topological properties - manifested by special surface states in the 
form of Fermi arcs, that connect the bulk Dirac points. This Dirac phase, a three dimensional analog 
of graphene, is proposed as the ground state of Y2lr207 and related compounds. A narrow window 
of magnetic 'axion' insulator, with axion parameter = ir, may also be present at intermediate U. 
An applied magnetic field induces ferromagnetic order and a metallic ground state. 



Previously, some of the most striking phenomena in 
soUds, such as high temperature superconductivity Ij and 
colossal magnetoresistance[2] were found in transition 
metal systems involving 3d orbitals, with strong electron 
correlations. Now it has been realized that in Ad and the 
5d systems, whose orbitals are spatially more extended, 
a regime of intermediate correlation appears. More- 
over, they display significant spin-orbit coupling, which 
modifies their electronic structure as recently verified in 
Sr2lr04 [3j. This is a largely unexplored domain, but 
already tantalizing new phenomena have been glimpsed. 
For example, in the 5d magnetic insulator, Na2lr03U, a 
disordered ground state persists down to the lowest mea- 
sured temperatures, making it a prime candidate for a 
quantum spin liquid[5]. 

It is known that strong spin-orbit interactions can lead 
to a novel phase of matter, the topological insulator :6]. 
However, the experimental candidates uncovered so far 
have weak electron correlations. Recently, it was realized 
that the iridates are promising candidates to realize topo- 
logical insulators [7], and that Iridium based pyrochlores 
in particular [8^, provide a unique opportunity to study 
the interplay of Coulomb interactions, spin-orbit coupling 
and the band topology of solids. 

The pyrochlore iridates, with general formula ^2^207, 
where A = Yttrium, or a Lanthanide element, will be 
the main focus of this work. Both the A and Ir atoms 
are located on a network of corner sharing tetrahedra 
[ni[TD]. Pioneering experiments [TT] on the pyrochlore iri- 
dates, revealed an evolution of ground state properties 
with increasing radius of the A ion, which is believed 
to tune electron correlations. While A =Pr is metallic, 
A =Y is an insulator as low temperatures. Subsequently, 
it was shown that the insulating ground states evolve 
from a high temperature metallic phase, via a magnetic 



transition [T2I [E]. The magnetism was shown to arise 
from the Ir sites, since it also occurs in A =Y, Lu, where 
the A sites are non-magnetic. While its precise nature 
remains unknown, the absence of a net moment rules out 
f erromagnetism . 

We show that electronic structure calculations can nat- 
urally account for this evolution and points to a novel 
ground state whose properties are described here. First, 
we find that magnetic moments order on the Ir sites in 
a non-colinear pattern with moment on a tetrahedron 
pointing all-in or all-out from the center. This structure 
retains inversion symmetry, a fact that greatly aids the 
electronic structure analysis. While the magnetic pat- 
tern remains fixed, the electronic properties evolve with 
correlation strength. For weak correlations, or in the 
absence of magnetic order, a metal is obtained, in con- 
trast to the interesting topological insulator scenario of 
Ref. [5] . With strong correlations we find a Mott insula- 
tor, with all-in/all-out magnetic order. However, for the 
case of intermediate correlations, relevant to Y2lr207, 
the electronic ground state is found to be an unusual 
Dirac semi-metal, with linearly dispersing Dirac nodes 
at the chemical potential. Indeed, this dispersion is anal- 
ogous to graphene|14|. but occurs inside a three dimen- 
sional magnetic solid. The small density of states leads to 
a vanishing conductivity at low temperatures. The Dirac 
fermions here are rather different from those in three di- 
mensional semi-metals such as elemental Bismuth, which 
are inversion symmetric and non-magnetic. Here, the 
Dirac fermions at a particular momentum are described 
by a handedness (which is left or right handed), and a 
two component wavefunction. They cannot be gapped 
unless they mix with a fermion of opposite handedness. 
In contrast, Dirac fermions in Bismuth have four com- 
ponent wavef unctions, no particular handedness, and are 



typically gapped. Such a three dimensional electronic 
structure has, to our knowledge, not been discussed be- 
fore. 

A key property of this Dirac semi-metal phase of two 
component Dirac fermions, is unusual band topology, 
reminiscent of topological insulators. Since the bulk 
fermi surface only consists of a set of momentum points, 
surface states can be defined for nearly every surface mo- 
mentum, and are always found to occur on certain sur- 
faces. They take the shape of 'Fermi arcs' in the surface 
Brillouin zone, that stretch between Dirac points. Hence 
we term this phase topological Dirac metal. 

We also mention the possibility of an exotic insulat- 
ing phase emerging when the Dirac points annihilate in 
pairs, as the correlation are reduced. This phase shows 
a topological magnetoelectric effect [T5]. captured by the 
magneto-electric parameter 9 = it, whose value is pro- 
tected by the inversion symmetry. Since it is analogous to 
the axion vacuum in particle physics THJ , so we call it the 
= n Axion insulator. Although our LSDA-I-U-I-SO cal- 
culations find that a metallic phase intervenes before this 
possibility is realized, we note that LDA systematically 
underestimates gaps, so this scenario could well occur in 
reality. Finally, we mention that modest magnetic fields 
could induce a reorientation of the magnetic moments, 
leading to a metallic phase. Our results are summarized 
in the phase diagram Figure [TJ Previous studies consid- 
ered ferromagnetism [T7], and structural distortion [TS] 
in iridates. 

Our calculations suggest that new functionalities such 
as controlling electrical properties via magnetic textures, 
and field induced metallic states can be realized in these 
materials, with implications for spintronics, magneto- 
electrical and magneto-optical devices. 
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FIG. 1: Sketch of the phase diagram for pyrochlore iridiates 
from our microscopic electronic structure calculation: Hori- 
zontal axis corresponds to the increasing interaction among 
Ir 5d electrons (the scale is obtained using LSDA-I-SO-I-U 
method) while the vertical axis corresponds to external mag- 
netic field which aligns the moments and triggers a transi- 
tion out of the zero field non-coUinear "all-in/all-out" ground 
state. We find normal magnetic metal, Mott insulator and 
Dirac semi-metal phases, and also possibly an exotic insu- 
lator, a magnetic insulator with magneto-electric parameter 
— TT, which we label here as axion insulator. 



Throughout, we exploit inversion symmetry which con- 
strains the phase diagram, by tracking wavefunction par- 
ities at time reversal invariant momenta. Near electronic 
phase transitions, a low energy k.p theory is developed 
to understand qualitative features of the neighboring 
phases. Finally, topological band theory based on mo- 
mentum space Berry connections is utilized in deducing 
the physical properties of the phases. 



METHOD 

We perform our electronic structure calculations based 
on local spin density approximation (LSDA) to density 
functional theory (DFT) with the full-potential, all- 
electron, linear-muffin-tin-orbital (LMTO) method [l9]. 
We use LSDA-fU scheme PO] to take into account the 
electron-electron interaction between Ir 5d electrons and 
vary parameter U between and 3 eV for Ir 5d electrons 
to see what effects the on site Coulomb repulsion would 
bring to the electronic structure of Iridates. In general, 
we expect that U can be somewhere between 1 and 2 
eV for the extended 5d states. When the A site is a 
rare earth element, we also add the Coulomb interac- 
tion for the localized 4/ electrons and use U = 6 eV. 
We use a 24x24x24 k-mesh to perform Brillouin zone 
integration, and switch off symmetry operations in order 
to minimize possible numerical errors in studies of vari- 
ous (non-)collinear configurations. We use experimental 
lattice parameters 12 in all set ups. 



MAGNETIC CONFIGURATION 

We first study magnetic configuration and discuss our 
results for ¥2^207. Since the strength of the spin or- 
bit (SO) coupling is large for Ir 5d electrons, and leads 
to insulating behavior in Sr2lr04|3], we perform the 
LSDA-|-U-fSO calculations. There are four Ir atoms in- 
side the unit cell forming a tetrahedral network as shown 
in Fig.lwhich is geometrically frustrated. Thus, we carry 
out several calculations with the initial state to be (i) 
ferromagnetic, with moment along (100), (HI), (HO) or 
(120) directions (ii) antiferromagnetic with two sites in 
a tetrahedron along and other two pointed oppositely to 
the directions above; non-colinear structures (iii) "all- 
in/out" pattern (where all moments point to or away 
from the centers of the tetrahedron see Fig 1), (iv) "2- 
in/2-out" (two moments in a tetrahedron point to the 
center of this tetrahedron, while the other two moments 
point away from the center, i.e. the spin-ice |21j configu- 




TABLE L The spin (S) and orbital (O) moment (in /is), 
the total energy Efot (in meV) for several selected magnetic 
configurations of Y2lr207 as calculated using LSDA+U+SO 
method with U=f .5 eV. The IDM is a coplanar configuration 
predicted for one sign of D-M interactions in Ref. [22j 

Configuration: (001) (111) 2-in/2-out IDM all-in/out 



FIG. 2: The pyrochlore crystal structure showing Ir corner 
sharing tetrahedral network and the magnetic configuration 
corresponding to the "all-in/all-out" alignment of moments. 
A degenerate state is obtained on reversing the moments. 



ration), and (v) "3~in/l-out" magnetic structures. 

We find that the "all-in/out" configuration is the 
ground state. Different from other magnetic configura- 
tions, during the self-consistency the "all-in/out" state 
will retain their initial input direction; thus, there is no 
net magnetic moment. This is consistent with the ab- 
sence of the magnetic hysteresis in experiments [T^. The 
all-in/all-out magnetic configuration was predicted to oc- 
cur in pyrochlore antiferromagnets with Dzyaloshinsky- 
Moriya (D-M) inter actions [22]. Symmetry dictates the 
form of D-M interactions except for the sign, which leads 
to two cases, direct and indirect D-M. The all-in/all-out 
state is the unique ground state for the former while the 
indirect D-M ground state is a coplanar state with the 
four spins being either antiparallel of orthogonal[22]- We 
find that the indirect D-M pattern has higher energy than 
the all-in/all-out state. This provides further evidence 
that the all-in/all-out spin configuration is the natural 
ground state. 

The next lowest energy configuration is the ferromag- 
netic state. Interestingly, the rotation of magnetization 
does not cost much energy despite strong SO interac- 
tions. The (111) direction is found to be lowest ferro- 
magnetic state, but the energy difference between this 
and the highest energy (001) state is just about 4.17 
meV per unit cell. Also, all of them produce a consider- 
able net magnetic moment in contrast to the experiment 
[TTl [T^ [53] . Our findings are summarized in Table I for a 
typical value of U=1.5 eV, and similar results are found 
for other values of U in the range from to 3 eV. We 
find that the energy difference between the ground and 
several selected excited states with different orientations 
of moments is small. Therefore, modest magnetic fields 
may induce a transition into the ferromagnetic state. 
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ELECTRONIC PHASES AND WAVEFUNCTION 
PARITIES 



We now discuss electronic properties of Iridates that 
emerge from our LSDA+U+SO calculations. A variety 
of phases ranging from normal metal at small U to Dirac 
semi- metallic at intermediate C/ ^ 1.5 eV and Mott in- 
sulating phase at U above 2 eV with non-coUinear mag- 
netic "all-in/out" ordering are predicted. Since pressure 
or chemical substitution may alter the screening and the 
electronic bandwidth resulting in changes in U we ex- 
pect that these phases can be observed experimentally in 
iridates. 

The basic features of the electronic structure can be 
understood by considering each of four Ir atoms in py- 
rochlore lattice which is octahedrally coordinated by six 
O atoms. This makes the Ir 5d state split into doubly de- 
generate Bg and triply degenerate t2g states. Due to the 
extended nature of Ir 5d orbital, the crystal-field split- 
ting between t2g and eg is large with the e^ band to be 
2 eV higher than the Fermi level. The bands near the 
Fermi level are mainly contributed by Ir t2g with some 
mixing with O 2p states. SO coupling has a considerable 
effect on these t2g states: it lifts their degeneracy and 
produces quadruplet with J^ff — 3/2 and doublet with 
Jeff = 1/2 with the latter lying higher in energy. Al- 
ternatively, one can view the d-shell in strong SO limit 
which splits the 10 fold degenerate state onto lower ly- 
ing J = 3/2 and higher lying J = 5/2 multiplets. A 
subsequent application of the cubic crystal field would 
leave the J = 3/2 multiplet degenerate but would split 
the 6-fold degenerate J = 5/2 state onto Fy doublet and 
Fg quadruplets. Since Ir occurs in its 4+ valence, its 5 
electrons would fill completely J — 3/2 subshell and put 
an additional one electron into the Fy doublet thought as 
the state with Je// = 1/2 in some previous work ^3,. For 
the solid, we thus expect 8 narrow energy bands at half- 
filling to appear in the vicinity of the chemical potential 
that correspond to the Fy doublets of four Ir atoms. 

The precise behavior of these electronic states depends 
on magnetic configuration. Our band structure calcula- 
tions for coUinear alignments of moments show metallic 
bands regardless the value of U that we use in our sim- 
ulations. On the other hand, we find that the electronic 
states for the non-coUinear "all-in/out" magnetic state 




TABLE 11: Calculated parities of states at Time Reversal In- 
variant Momenta (TRIMs) for several electronic phases of the 
iridates. Only the top four filled levels are shown in order of 
increasing energy. 



FIG. 3: Evolution of electronic band structure of Y2Ir207 
shown along high symmetry directions, calculated using 
LSDA+U+SO method with three different values of U equal 
(a) 0, (metallic) (b) 1.5 eV, and (c) 2 eV. (Insulator with 
small gap). The Dirac point that is present in case (b), is not 
visible along high symmetry lines. 



depend strongly on the actual value of U used in the cal- 
culation. In particular, we predict that when U is less 
than 1 eV, the ground state is a normal metal while if 
U is about 1.8 eV or larger, we find the band structure 
to be insulating with an energy gap whose value depends 
on U. 




FIG. 4: Calculated energy bands in the vicinity of the Fermi 
level using LSDA-|-U-fSO method with U=1.5 eV. on left: 
corresponds to plane kz — with band parities shown; (b) 
corresponds to plane kz — 0.32n/a where Dirac point is pre- 
dicted to exist. The shaded plane is at the Fermi level. 

Weak Correlations: An interesting recent study pro- 
posed a tight-binding model for the non-magnetic phase 
of the iridates, which was a topological insulator [H], a 
natural phase on the pyrochlore lattice [HI [H] • Our LDA 
studies of the realistic electronic structures contradict 
this, instead we find a metallic phase (see Fig. Isk). One 
can understand the discrepancy by analyzing the struc- 
ture of energy levels at the F point (Brillouin Zone center) 
for the low energy 8-band complex, composed of the four 
Jeff = 1/2 states. In ^, these appear with degeneracies 
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4,2,2 which after filling with 4 electrons results in an in- 
sulating band structure. Our study of the non-magnetic 
state using LDA-t-SO method (with no U) results, on 
the other hand, in the sequence 2,4,2 of degeneracies, 
which is necessarily metallic assuming 4 levels are filled. 
Thus, one needs either magnetic order, or a structural 
transition^lSj. to recover an insulating phase. 

Strong Correlations and the Mott Limit: When U > 
l.SeV, an insulating band structure is obtained with the 
all-in/all-out magnetic configuration, as shown in the fig- 
ure |3t. Indeed, it remains qualitatively similar on in- 
creasing U to large values, where a site localized mo- 
ment is expected, i.e. the Mott insulator. This can 
be further verified by calculating the parity eigenvalues. 
Note that all the magnetic structures considered above 
preserve inversion (or parity) symmetry. Thus, if we 
pick e.g. an Iridium atom as the origin of our coor- 
dinate system, then inversion r — >■ — r leaves the crys- 
tal structure and magnetic pattern invariant. This im- 
plies a relation between crystal momenta ±k. At special 
momenta, called TRIMs (Time Reversal Invariant Mo- 
menta), that are invariant under inversion, we can label 
states by parity eigenvalues ^ = ±1. In the Brillouin zone 
of the FCC lattice these correspond to the F = (0, 0, 0), 
and X, Y, Z [=27r(l, 0, 0) and permutations] and four 
L points [7r(l, 1, 1) and equivalent points]. These par- 
ities are very useful to study the evolution of the band 
structure. The TRIM parities of the top four occupied 
bands, in order of increasing energy, are shown in Table 
2. Note, although by symmetry all L points are equiv- 
alent, the choice of inversion center at an Iridium site 
singles out one of them, L' . With that choice the pari- 
ties at the two sets of L points are the opposite of one 
another. The parities of the all-in/out state remains un- 
changed above U > Uc ^ 1-8 eV, and is shown in the top 
row under U = 2 eV. It is readily seen that these parities 
arise also for a site localized picture of this phase, where 
each site has an electron with a fixed moment along the 
ordering direction. Due to the possibility of such a lo- 
cal description of this magnetic insulator, we term it the 
Mott phase. 

Intermediate Correlations: For the same all-in/out 
magnetic configuration, at smaller U — 1.5 eV, the band 
structure along high symmetry lines also appears insulat- 
inglsb, and at first sight one may conclude that this is an 
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FIG. 5: Brillouin zone of the FCC lattice with locations of 
Dirac points (shown by + signs denoting their "positive" chi- 
ral charges) as found by our LSDA+U+SO calculation with 
U=1.5 eV for Y2Ir207. 



extension of the Mott insulator. However, a closer look at 
the parities reveals that a phase transition has occurred. 
A pair of levels with opposite parity are exchanged at 
the L points. Near this crossing point it can readily be 
argued that only one of the two adjacent phases can be 
insulating[30]. Since the large U phase is found to be 
smoothly connected to a gapped Mott phase, it is reason- 
able to assume the smaller U phase is the non-insulating 
one. This is also borne out by the LSDA+U+SO band 
structure. A detailed analysis perturbing about this tran- 
sition point (the k.p expansion see [30]) allows us to show 
that a Dirac semi-metal is expected for intermediate f/, 
with 6 Dirac nodes about every L point. Indeed, in the 
LSDA+U+SO band structure, we find a 3 dimensional 
Dirac crossing located within the FXL plane of the Bril- 
louin zone. This is illustrated in Fig[5] and corresponds 
to the k-vector (0.52,0.52, 0.3)27r/a. There also are two 
additional Dirac points in the proximity of the point L 
related by symmetry. When U increases, these points 
move toward each other and annihilate all together at 
the L point close to t/ = 1.8 eV. This is how the Mott 
phase is born from the Dirac phase. Since we expect that 
for Ir 5d states the actual value of the Coulomb repulsion 
should be somewhere within the range leV < U < 2eV 
we thus conclude that the ground state of the Y2lr207 is 
most likely the semi-metallic state with the Fermi surface 
characterized by a set of Dirac points but in proximity 
to a Mott insulating state. Both phases can be switched 
to a normal metal if Ir moments are collinearly ordered 
by a magnetic field. 



Topological Dirac semi-metal 

The effective Hamiltonian in the vicinity of a Dirac 
momentum k = ko + q is: 



Hr 



^ v., • qcTi 



(1) 



, where energy is measured from the chemical poten- 
tial, and ai are the three Pauli matrices. The three ve- 
locity vectors v^ are generically non-vanishing and lin- 
early independent. The energy dispersion is AE = 

^\/X]j=i(^j • q)^. Note, by inversion symmetry, there 
must exist Dirac points at both ko and —kg, whose ve- 
locity vectors are reversed. One can assign a chirality 
(or chiral charge) c = ±1 to the fermions defined as 
c — sign(vi • V2 X V3), so Dirac points related by in- 
version have opposite chirality. Note, since the 2x2 
Pauli matrices appear, our Dirac particles are two com- 
ponent fermions. In contrast to regular four component 
Dirac fermions, it is not possible to introduce a mass gap. 
The only way to eliminate them is if they meet with an- 
other 2 component Dirac dispersion in the Brillouin Zone, 
but with opposite chiral charge. Thus they are topo- 
logical objects. We note that near each L point there 
are three Dirac points related by the three fold rotation, 
which have the same chiral charge. Fig[5] denotes those 
points as "+" dots. Another three Dirac points with op- 
posite chirality, related by inversion. Thus, there arc 24 
Dirac points in the whole Brillouin zone. Since they are 
all related by symmetry, they are at the same energy. 
The chemical potential is fixed to be at the Dirac point 
energy as verified in the microscopic calculation. The 
Fermi velocities at the Dirac point are found to be typ- 
ically an order of magnitude smaller than in graphene. 
We briefly note that this Dirac semimetal is a critical 
state with power law forms for various properties, which 
will be described in more detail elsewhere. For example, 
the density of states N{E) ex E'^. The small density of 
states makes it an electrical insulator at zero tempera- 
ture. For a single node with isotropic velocity v, the a.c. 
conductivity in the free particle limit of the clean system 
is ain) = ^^-ftanh\n\/ AksT. 

Surface States: We now discuss surface states that 
are associated with the presence of the two component 
Dirac fermions. We first note that they behave like 
'magnetic' monopoles of the Berry fiux whose charge 
is given by the chirality.The Berry connection, a vec- 
tor potential in momentum space, is defined by .A(k) = 
J2n=i *(wnk| Vk|u„k) where N is the number of occupied 
bands. As usual, an analog of the magnetic field, the 
Berry flux, is defined as J^ = Vk x A. Now consider 
energy eigenstates at the Fermi energy (taken to be at 
E = 0). In the bulk, this corresponds to the set of Dirac 
points, hence the bulk Fermi surface is a collection of 
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FIG. 6; Illustration of surface states arising from bulk Dirac 
points. For simplicity, only a pair of Dirac points with oppo- 
site chirality are shown. The imaginary cylinder in momen- 
tum space has unit Chern number, due to the Berry monopole 
at the Dirac point. Hence a surface state must arise, as shown 
schematically in the same plot. When the Fermi energy is at 
the Dirac point, a Fermi arc is present which connects the sur- 
face momenta of the projected bulk Dirac points of opposite 
chirality. 



Fermi points. However, in the presence of a surface (say 
the plane z = 0), new low energy states may be gener- 
ated. We show that these will occur along a curve in the 
surface Brillouin zone as is illustrated in Fig. 5. The end 
points of this curve occur at the bulk Fermi point mo- 
menta, projected onto the surface Brillouin Zone. Also, 
the curve connects Dirac nodes with opposite monopole 
charge. If more than one Dirac node projects to the same 
surface momentum, the sum of the monopole charges 
should be considered. This is argued by showing that 
there must be Fermi arcs on the surface Brillouin zone 
emanating from the projection (fcox, koy) of the monopole 
as argued below. 

Origin of Surface States: We now prove that the band 
topology associated with the Dirac point leads to sur- 
face states. Construct a curve in the surface Brillouin 
zone encircling the projection of the bulk Dirac point, 
which is traversed counterclockwise as we vary the pa- 
rameter A : — ?► 27r; k^ — {kx{\),ky{\)) (see Fig. |6]). 
We show that the energy e\ of a surface state at momen- 
tum kx crosses E = 0. Consider H{X,kz) = H{kx,kz), 
the gapped Hamiltonian of the two dimensional subsys- 
tem defined by this curve. The two periodic parameters 
A, kz define the surface of a torus in momentum space 
(see Fig. [6]) . The Chern number of this two dimensional 
band structure is given by the Berry curvature integra- 
tion: ^ / TdkzdX which, by Stokes theorem, simply cor- 
responds to the net monopole density enclosed within 
the torus. This is obtained by summing the chiralities 
of the enclosed Dirac nodes. Consider the case when the 



net chirality is unity, corresponding to a single enclosed 
Dirac node. Then, the two dimensional subsystem de- 
fines a quantum Hall insulator with unit Chern number. 
When defined on the half space z < 0, this corresponds to 
putting the quantum Hall state on a cylinder, and hence 
we expect a chiral edge state. Its energy e^ spans the 
band gap of the subsystem, as A is varied. Hence, this 
surface state crosses zero energy somewhere on the sur- 
face Brillouin zone k^o . Such a state can be obtained for 
every curve enclosing the Dirac point. Thus, at zero en- 
ergy, there is a Fermi line in the surface Brillouin zone, 
that terminates at the Dirac point momenta (see Fig. 
p|. An arc beginning on a Dirac point of chirality c has 
to terminate on a Dirac point of the opposite chirality. 
Clearly, the net chirality of the Dirac points within the 
(A, kz) torus was a key input in determining the number 
of these states. If Dirac points of opposite chirality line 
up along the kz direction, then there is a cancelation and 
no surface states are expected. 

For U = 1.5 eV, a Dirac node is found to occur 
at the momentum (0.52, 0.52, 0.31)27r/a and equivalent 
points (see Fig. 3). They can be thought of as occur- 
ing on the edges of a cube, with a pair of Dirac nodes 
of opposite chirality occupying each edge, as, e.g., the 
points (0.52, 0.52, 0.31)27r/a and (0.52,0.52, -0.31)27r/a. 
For the case of C/ — 1.5 eV, the sides of this cube have 
the length 0.52(47r/a). Thus, the (111) and (110) sur- 
faces would have surface states connecting the projected 
Dirac points. If, on the other hand we consider the sur- 
face orthogonal to the (001) direction, it would lead to 
the Dirac points of opposite chirality being projected to 
the same surface momentum, along the edges of the cube. 
Thus, no protected states are expected for this surface. 

Model Calculation: To verify these theoretical con- 
siderations, we have constructed a tight binding model 
which has features seen in our electronic structure cal- 
culations for YIr207. We consider only t2g orbitals of Ir 
atoms in the global coordinate system. Since Ir atoms 
form tetrahedral network (see Fig. [2]), each pair of 
nearest neighboring atoms forms a corresponding cr— like 
bond whose hopping integral is denoted as t and other 
two TT— like bonds whose hopping integrals are denoted 
as t' . To simulate the appearance of the Dirac point it 
is essential to include next-nearest neighbor interactions 
between t2g orbitals which are denoted as t" . With the 
parameters t = 0.2, t' = 0.5i, t" = — 0.2<, the value of the 
on-site spin-orbit coupling equal to 2.5i and the applied 
on-site splitting between spin up and spin down states 
equal to 0.1 referred to the local quantization axis which 
simulates our non-coUinear 'all-in/out' configuration we 
can model both the bulk Dirac metal state and its sur- 
face. The calculated (110) surface band structure for the 
slab of 128 atoms together with the sketch of the obtained 
Fermi arcs is shown in Fig. [7j Notice that since the slab 
calculation involves two surfaces, the corresponding sur- 
face states and Fermi arcs for both surfaces are generated. 
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FIG. 7: Calculated surface energy bands corresponding to 
(110) surface of the pyrochlore iridate Y2Ir207. A tight bind- 
ing approximation has been used to simulate the bulk band 
structure with 3D Dirac point as found by our LSDA+U+SO 
calculation. The plot corresponds to diagonalizing 128 atoms 
slab with two surfaces. The top inset shows the deduced Fermi 
arcs connecting projected bulk Dirac points of opposite chi- 
rality. The bottom inset shows a sketch how these Fermi arcs 
are expected to behave for the (111) surface. 



We also display the expected surface states for the (111) 
surface. Note, no special surface states are expected for 
the (001) surface. 



OTHER TOPOLOGICAL PHASES 

We recall that topological insulators are non-magnetic 
band insulators with protected surface states [6]. Time 
reversal symmetry is required in the bulk to define 
these phases. When the surface states are eliminated 
by adding, for example magnetic moments only on the 
surface, a quantized magneto-electric response is ob- 
tained, where a magnetic field induces a polarization: 
P — 9^^'B, with the coefficient \9\ [TS] is only defined 
modulo 2tt. Under time reversal, 6 — > —9. Apart from 
the trivial solution 9 = 0, the ambiguity in the defini- 
tion of 9 allows also for 9 = ir. For topological insulators 

9 = 77. 

In magnetic insulators, 9 is in general no longer quan- 
tized [26j . However, when inversion symmetry is re- 
tained, 9 is quantized again, since inversion also changes 
its sign. Thus again 9 = 0, 7rmod27r, and an insulator 
with the latter value will be termed a, 9 = ir axion insu- 
lator. 

Which is the appropriate description of the pyrochlore 



iridate phases we have described? For the Mott insulator, 
at large U, the charge physics must be trivial and so we 
must have 9 = 0. Next, since the Dirac semi metal phase 
is gapless in the bulk, 9 is ill defined. However, we note 
that on reducing U, the location of the Dirac points shift, 
with nodes of opposite chirality approaching each other. 
If these meet and annihilate, then one recovers a gapped 
phase in the low U regime. However, in the process the 
resulting phase will have 9 = n. Indeed, the presence of 
the intervening Dirac phase can be deduced from the re- 
quirement that 9 has to change between these two quan- 
tized values. As described elsewhere|25j. the condition 
for 9 = 7: when deduced from the parities, turns out to 
be the same as the Fu-Kane formula, for time reversal 
symmetric insulators [571 121] > i-e. if the total number 
of filled states of negative parity at all TRIMs taken to- 
gether is twice an odd integer, then 9 = tt. Otherwise 
9 = 0. The small U insulator has the same parities as 
the Dirac semi-metal, since the Dirac points annihilate 
away from a TRIM. From Table [III we can see that indeed 
this corresponds to 9 = tt, since there are 14 negative par- 
ity filled states, while the Mott insulator corresponds to 
9 = 0, having 12 negative parity filled states. 

Unfortunately, within our LSDA+U+SO calculation, a 
metallic phase intervenes on lowering U < l.OeV, before 
the Dirac points annihilate to give the axion insulator. 
We point out this possibility nevertheless, since LDA sys- 
tematically underestimates the stability of such gapped 
phases. Moreover it provides an interesting example of 
a pair of insulators with the same magnetic order, but 
which are nevertheless different phases. Thus, the Mott 
insulator is distinct from the smaller U 'Slater' insulator, 
unlike in many other cases where they are smoothly con- 
nected to one another. Inversion symmetry is critical in 
preserving this distinction. 

In summary, a theoretical phase diagram for the phys- 
ical system is shown in Figure [l| as a function of U and 
applied magnetic field, which leads to a metallic state 
beyond a critical field. The precise nature of these phase 
transformations are not addressed in the present study 
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APPENDIX 



I. Effective k.p Theory and Intervening Dirac Metal 
Phase 



induce a matrix element between the levels. And 
also, since kj^ is a 2 vector, it transforms under the 
rotation - the rotationally symmetric form allowed 
is AH = k^ cos36'tj. + k^ sm36Ty. Putting this all 
together we have the effective Dirac Hamiltonian 
near the L point: 



Consider a pair of states at the L point which have 
opposite parity, and cross each other as we tune U. We 
want to understand what happens to the band structure. 

The L point has two symmetries which do not change 
its crystal momentum. First of course is inversion, and we 
can label states by the eigenvalues P — ±1. The second is 
120° rotations about a line joining L — T. There are three 
possible eigenvalues which we call s — —1/2, 1/2, 3/2. 
So, any state at this point can be labeled by {P, s}. Now 
consider writing the effective 2x2 Hamiltonians for a pair 
of states that are near each other in energy: 

1. At the L point: Since we have inversion symmetry, 
the two states taken to be eigenvalues of t^ — ±1, 
cannot mix. Hence the effective Hamiltonian is 

H{L) =. At, 

where the coefficient A changes sign when the lev- 
els pass through each other. Note, the s quantum 
number of the two levels is irrelevant here. 

2. Along the F — L Direction: We still have the quan- 
tum number s, but not P, since that inverts the 
momentum. Denoting by fc, the momentum along 
this line deviating from the L point, we have two 
cases. If the s quantum number of the two levels 
is different, they still cannot mix, so the effective 
Hamiltonian is iJ = (A -I- k1)Tz- Now, when A < 
there are two solutions fc, = ±\A^A, where there 
are nodes along this T — L line. You can see this 
for the s = 1/2 and s — 3/2 crossings in the data. 

However, if the states have the same s quantum 
number they can mix, once you move away from L. 
Now the effective Hamiltonian is: 



H{T -L) = {A 



/c?- 



where the second term arises since inversion is bro- 
ken on moving away from L allowing for mixing. 
Now, the spectrum is i<^ = ^^(fc^ — |A|)2 -|- k^, for 
A < 0, so despite a level crossing there is no node 
along the T — L line. 

General Point in BZ: In the latter case, does this 
mean there are no Dirac points? No - we just need 
to move off the T — L line. Let the deviation be 
k_L, a 2 vector. The fact the 2 levels have opposite 
parity means we need an odd function of k±^ to 



iJ(k) = (A + fc^)r, + (fc, + k\ cqsW)t,^ + k\ smWTy 

Note, this has the form A{k)Tz + B(k)Tx + C{k)Ty. 
For a node, A = B = C = Q. This occurs \i:A = 
so fcz = +-\/— |to|, and C = implies 9 = P^^/'i^ 
where p is an integer; and finally k\ — kz, when 
p — 1,3,5. Similarly for kz < 0, the nodes are 
inverted. In all we have 6 nodes for this L point, 
24 Nodes in all. 

Actually this is not the complete expansion. Strictly 
we should write A{k) = m + k^ + ak\. This is an ef- 
fective mass Hamiltonian near the L point. If a > 0, 
then it turns out the conclusions are the same as above, 
without this additional term (i.e., there is a transition 
from 6 = 7: magnetic axion insulator to Dirac metal on 
increasing U). However, if a < 0, this completely changes 
the conclusions, as discussed below. In fact, this turns 
out to be the physically relevant case according to the 
electronic structure calculations for Y2lr207. 

Let us assume a = —l/m2 < 0. Then the effective 
Hamiltonian: 



H{k) = (A- 



fci 



-{l3kz+kl^ cos39)Tx+kl^ sin36'Ty 



2?Til 2?Ti2 ' 

where a few parameters have been labeled. The Dirac 
nodes then are at: (i) C = — > 6* = pir/S and (ii) using 
this, i? = — > fcz = ztkj_/l3 depending on whether we 
look at p =odd or p = even. Finally, using these relations 
we have for A = equation: 



A + 



kl 



kl 



2l3^mi 2?Ti2 



= 



for small A this has the solution k'j_ = 2m2A. Note, 
this has a solution only for A > 0, i.e. before the gap gets 
inverted at the L point on increasing U. Thus, in this 
scenario, there is a Dirac point only in the small U phase. 
These Dirac points live along = pw/3 plane, which is 
the kx = ky plane that contains the points T — L — K 
and rotations thereof. 



II. Summary of Experiments 

We now summarize the experimental facts about the 
pyrochlore iridates j42lr207. Early work revealed that 
increasing the A ion radius triggered a metal-insulator 



transition in the ground state. Thus while A =Pr is 
metallic, A =Y is insulating at low temperatures [11]. 
Reducing the ionic radius is believed to narrow the band- 
width and increase correlations. Subsequent improve- 
ment in sample quality revealed that A =Eu, Sm, Nd 
also displayed low temperature insulating states [13]. In 
these systems, a metal insulator transition is clearly ob- 
served on cooling (eg. at Tmi = 120K for A =Eu). At 
the same temperature, a signature in magnetic suscepti- 
bility is also observed, indicative of a magnetic transition 
[12j . This magnetic signature, wherein the field cooled 
and zero field cooled magnetic susceptibilities separate 
below the transition temperature, is reminiscent of a spin 
glass state. Since this signature is see in A =Y, Lu [12], 
with nonmagnetic A site atoms, it is associated with Ir 
site moments. In Y2lr207, no sharp resistivity signature 
has been reported at the magnetic transition, but the 
resistivity climbs steeply on cooling below this temper- 
ature. Moreover, light hole doping suppresses both the 
insulating state and the magnetic transition [23] . Finally, 
we note that a thermodynamic signature of the magnetic 
transition, a bump in the specific heat, is observed in 
clean samples of A =Sm [T^]. X-ray scattering did not 
observed any structural change below the ordering tran- 
sition [13] , although the presence of new lines in the Ra- 
man spectroscopy j^Sj has been attributed to lowering of 
cubic symmetry m A = Eu, Sm but not in A —Nd. 

We now discuss how the present theoretical description 
sits with these facts. We propose that the low tempera- 
ture state of Y2lr207 (and also possibly oi A = Eu, Sm 
and Nd iridates) is the Dirac semi-metal, with all-in/all- 
out magnetic order. This is broadly consistent with the 
interconnection between insulating behavior and mag- 
netism observed experimentally. It is also consistent with 
being proximate to a metallic phase on lowering the cor- 
relation strength, such as A =Pr. In the clean limit, a 
three dimensional Dirac semi-metal is an electrical insu- 
lator, and can potentially account for the observed elec- 
trical resistivity. The noncolinear magnetic order pro- 
posed has Ising symmetry and could undergo a contin- 
uous ordering transition. Since this configuration is not 
frustrated, it is not clear how spin glass behavior would 
arise, but the observed magnetic signature could perhaps 
also arise from defects like magnetic domain walls. A di- 
rect probe of magnetism is currently lacking and would 
shed light on this key question. 
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